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Mapping the shores of the brown dwarf desert. IV. Ophiuchus 


Anthony C. Cheetham^’^’*, Adam L. Krans^, Michael J. Ireland^, Lucas Cieza®’®, Aaron C. 

Rizzuto^, Peter G. Tuthilh 


ABSTRACT 

We conduct a multiplicity survey of members of the p Ophiuchus cloud complex with high 
resolution imaging to characterize the multiple star population of this nearby star forming region 
and investigate the relation between stellar multiplicity and star and planet formation. Our 
aperture masking survey reveals the presence of 5 new stellar companions beyond the reach of 
previous studies, but does not result in the detection of any new substellar companions. We find 
that 43±6% of the 114 stars in our survey have stellar mass companions between 1.3-780 AU, 
while 7^5% host brown dwarf companions in the same interval. By combining this information 
with knowledge of disk-hosting stars, we show that the presence of a close binary companion 
(separation < 40 AU) significantly influences the lifetime of protoplanetary disks, a phenomenon 
previously seen in older star forming regions. At the ~l-2Myr age of our Ophiuchus members 
~2/3 of close binary systems have lost their disks, compared to only ~30% of single stars and 
wide binaries. This has significant impact on the formation of giant planets, which are expected 
to require much longer than 1 Myr to form via core accretion and thus planets formed via this 
pathway should be rare in close binary systems. 

Subject headings: binaries: general — brown dwarfs — stars: low-mass — stars: pre-main sequence 


1. Introduction 

Multiplicity surveys provide some of the most 
stringest tests of stellar formation theories. The 
distribution and frequency of companions are key 
predictions of these theories that are easily tested. 
Radial velocity (RV) surveys have had great suc¬ 
cess at revealing close stellar and substellar com¬ 
panions, especially around older stars, while high 
resolution imaging surveys have increased our un- 
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derstanding of companions at wider separations. 

RV surveys have revealed a wealth of both stel¬ 
lar and planetary mass companions in close orbits 
around their host stars. However, they have re¬ 
vealed a surprising lack of brown dwarf mass com¬ 
panions in close orbits, with estimates suggest- 


ing frequencies of < 1% (1 

darcy & Butler 

2000 

Grether & Lineweaver 

200 

31. This phenomenon 


has been labelled the “brown dwarf desert”. 


In contrast, imaging surveys targeting wider 
separations have found that the frequency of 
such companions may not be anomalously low, 
but rather an extension of the binary mass-ratio 


function to lower masses (Metchev & Hillenbrand 


20091. 


In this paper we investigate the crucial sepa¬ 
rations between these two approaches, following 
on from previous work in Kraus et al. (2008) and 


Kraus et al. (2011). These studies investigated the 


Upper Scorpius subgroup of the Sco-Cen OB asso¬ 
ciation and the Taurus-Auriga star forming region 
respectively. Both studies resulted in companion 
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distributions consistent with a flat mass ratio dis¬ 
tribution, finding 5 and 6 companions with mass 
ratios q< 0.1 respectively. This provides further 
evidence that the frequency of brown dwarf com¬ 
panions is consistent with the observed trend for 
stellar mass companions. 


The distributions of orbital parameters in bi¬ 
nary stars offers another key test of star formation 
models, through comparison of their predictions to 
the observed properties of multiple systems. Sev¬ 
eral studies have targeted binaries in the field, 
finding evidence that these properties are mass- 


dependent ( Duquennoy &: Mayor| 1991 Fischer & 
Marcyl [T99^ Close et al. |2003| [Raghavan et al. 

20101. For example, companions to Solar mass 


stars have a higher mean separation than compan¬ 
ions to low mass stars. In addition, the number 
of low mass companions to Solar mass stars ap¬ 
pears high, while the distribution of companion 
mass ratios for low mass stars is peaked more to¬ 
wards unity. 

However, surveys of young star forming regions 
have shown significant differences in many of these 
properties. The overall binary frequencies in these 
regions appear to be much higher than the fleld 


(> 80% 

Chez et al. 1993 

Simon et al. 1995 

Kohler 

et al. 

2000 

Kraus et al. 

2008 

2011), and the 


shapes of the separation and mass ratio distri¬ 
butions appear different. Dynamical interactions 
may play a significant role in causing these dif¬ 
ferences, since studies of denser regions such as 
young clusters have shown similar results to the 
fleld ( Duchene et al.|[l999 Kohler et al.||2006 ). A 
large fraction of stars in the field are thought to 
originate from these dense clusters, and the simi¬ 
larities between their properties echoes this idea. 


Samples less affected by dynamical interactions 
provide simpler tests of binary formation pro¬ 
cesses, and so obtaining robust statistics for young 
star forming regions like Ophiuchus are important 
to test multiple star formation theories. 


The relationship between binarity and disk evo¬ 
lution is another key link to understanding star 
and planet formation. While a large body of work 
has concentrated on the evolution and formation 
of planets in single star systems, the majority of 
solar type stars exist in multiple systems which 
may have a profound effect on the way in which 
these processes occur. 


Surveys targeting disks and stellar multiplicity 
in nearby star forming regions have shown correla¬ 
tions between the presence of a binary companion 
and the properties and presence of a circumstellar 


disk 

Ghez et aL|1997t 

Cieza et al.||2009 

Duchene 

2010 

Kraus et al.|2012 

1. These studies found that 


the presence of a close (< 40 AU) binary compan¬ 
ion can significantly speed up the dispersal or in¬ 
hibit the formation of protoplanetary disks. 


By combining the results of previous multiplic¬ 


ity surveys of several star forming regions, Cieza 


et al. (20091 found that close binaries with separa¬ 


tions less than 40 AU were half as likely to retain 
their disks as binary systems with larger separa¬ 
tions. However, the timescale of this effect remains 
unclear. 


Comparing the results for different star form¬ 
ing regions by age, [Kraus et al. (20121 found that 
~2/3 of all close binaries have no disk at ages 
of 1-2 Myr in Taurus. Despite this, stable con¬ 
figurations appear to exist that allow some disks 
around close binaries to persist for ~ 10 Myr in 
Upper Scorpius. In contrast, the presence of a 
wider companion does not appear to affect the life¬ 
time of protoplanetary disks at these ages. Com¬ 
parison of results from the 1-2 Myr old Ophiuchus 
region with the regions investigated in previous 
studies will provide information on the age depen¬ 
dent properties of these effects and the scatter be¬ 
tween individual regions. 

The relationship between the timescales of disk 
dispersal and giant planet formation defines the 
relative abundance of such planets. The two 
canonical giant planet formation theories of core 
accretion ( Pollack et al.|1996 l and disk instability 
( Boss|[20^ I predict different timescales. Core ac¬ 
cretion requires several Myr to form giant planets 
from protoplanetary disks ( Hubickyj et al.|[2005 l, 
while disk instability is most efficient at much 
younger ages (<0.5Myr, Boss 2001). The short 
lifetime of disks in close binary systems would 
lead to a low occurence rate of giant planets that 
formed through core accretion, a key prediction to 
differentiate between the two theories. 


We utilise the technique of Sparse Aperture 
Masking (SAM) to perform a high-resolution sur¬ 
vey of the nearby Ophiuchus star forming region. 
The relatively small distance to this association as 
well as its relative youth provide an ideal opportu¬ 
nity to investigate changes in both the prevalence 
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of brown dwarf companions and the relationship 
between disk evolution and multiplicity. Compar¬ 
ison of results from the ~l-2Myr old Ophiuchus 
region with the older regions investigated in pre¬ 
vious studies will provide information on the age 
dependent properties of these effects. 

2. Survey Sample 



B5 AO A5 FO F5 GO G5 KO K5 MO M5 


Spectral Type 

Fig. 2.— The distribution of spectral types for the 
input sample described in Section]^ 


This survey consists of observations of mem¬ 
bers of the p Ophiuchus cloud complex. The star 
formation history of this region is complex, with 
several periods of recent star formation. Its similar 
location and distance to the Sco-Cen OB associa¬ 
tion has led to the suggestion that the most recent 
episode of star formation in this region was trig¬ 
gered by an interaction between the large L1688 
cloud and a shock wave from the Upper Scorpius 
sub group, occuring approximately 1-1.5 Myr ago 
( Vrba|[l9^ Loren fc Wootten||1986 1. 

The first estimates of the distance to the 
cloud complex suggested a value of around 160 pc 
(Bertiau 1958 Whittet 1974). However, more 
recent studies have measured the distance to be 


between 119-135pc (de Geus et al. 1989 Lom¬ 


bardi et al. 2008 Mamajek 2008). In this study. 


we adopt a value of 130 pc for all of our association 
members. 


Members of Ophiuchus have a wide range of 
ages consistent with several episodes of star for¬ 
mation in different regions. A young population of 
Young Stellar Objects (YSOs) in the centre of the 


cloud with a median age of 0.3 Myr was identified 


by Greene & Meyer (19951 and Luhman & Rieke 


(1999), but an older population of YSOs cover¬ 
ing a larger area has an age of 2 Myr, as noted by 


Wilking et al. (20051. These studies estimated the 


age of observed members through the comparison 
of H-R diagram positions to model evolutionary 
tracks, a process that can introduce large system¬ 
atic effects, as illustrated by the analysis of F-type 
stars in Sco-Cen by Pecaut et al. (2012). In this 
study we adopt an age of 1-2 Myr and focus on 
members close to the L1688 cloud core. 


Our target list was arrived at by considering 
membership of the p Ophiuchus cloud complex. 
Since no canonical survey exists, a conservative 
approach was taken where all possible members 
from the literature were considered and then sub¬ 
ject to a cut in right ascension and declination 
(16 20 00 to 16 40 00 and -23 00 00 to -27 00 00 
respectively). This was chosen to cover a wide re¬ 
gion around the L1688 cloud that hosts ongoing 
star formation. The initial target list was made 


by compiling targets from Herbig & Bell 

19881, 

Bouvier & Appenzeller (1992), Meyer et al. 

19931, 

Martin et al. 

(1998), 

Preibisch et al. 

(1998 

, Luh- 

man & Rieke 

(1999) 

Wilking et al. 

2005) 

Cieza 

et al. (2010), 

and Erickson et al. (2011). 


Our targets were further restricted to those 
members with measured photometry and known 
spectral types. Considering the limiting magni¬ 
tudes for reasonable AO correction at Keck and 
the VLT, targets fainter than R = 15 mags and K 
= 9.5 mags (for the Keck NIRC2 visible WFS and 
the VLT NACO IR WFS respectively) were not in¬ 
cluded. In an attempt to give a sample with better 
mass completeness, stars with spectral types later 
than M4 were observed but excluded from statis¬ 
tical analysis. 

Out of the 236 possible Ophiuchus members 
from the literature, 114 fit these criteria. These 
targets are listed in Table and their distribu¬ 
tion on the sky is shown overlaid on an extinction 
map of the Ophiuchus L1688 cloud in Figure 
The distribution of spectral types for these tar¬ 
gets is shown in Figure This target list com¬ 
prises the input sample for our study, and consists 
of observed targets, known binaries that were not 
reobserved, and targets that were not observed. 

We have compiled the results of previous multi¬ 
plicity surveys of p Ophiuchus members performed 
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Fig. 1.— Distribution of targets that fit the criteria described in Section plotted over contours marking 
1, 2.5, 5, 10 and 15 mag of (B-V) extinction from Schlegel et al. ( 1998[ ). Several prominent nearby stars 
are marked with triangles and labelled for reference, as well as HD 147889, the highest mass target in our 
sample. Targets marked with red circles were not obseved with aperture masking. 


Table 1 
Survey Sample 


Name 

RA 

(J2000) 

Dec 

(J2000) 

SpT 

Mass 

{M@) 

R 

(mag) 

K 

(mag) 

Reference ^ 

Disk? 

Multiple 

System? 

DoAr 13 

16 20 39.60 

-26 34 28.4 

M2 

0.61 

14.89 

9.921 

6 

Y 

N 

GSC 6794-480 

16 20 45.96 

-23 48 20.8 

K3 

2.07 

12.66 

8.927 

7 

N 

N 

[MMG98] RX ,11620.9-2352 

16 20 57.87 

-23 52 34.3 

K3 

2.07 

9.88 

8.393 

6 

N 

N 

[MMG98] RX J1621.2-2342a 

16 21 14.50 

-23 42 20.0 

K7 

0.98 

12.51 

8.99 

6 

N 

N 

HIP 80126 

16 21 19.18 

-23 42 28.7 

B5 

5.90 

6.98 

6.418 

12 

Y 

N 

[MMG98] RX J1621.4-2312 

16 21 28.44 

-23 12 11.0 

K7 

0.98 

13.55 

8.739 

6 

N 

N 

[MMG98] RX J1621.4-2332b 

16 21 28.81 

-23 32 38.9 

MO 

0.82 

11.02 

7.167 

6 

N 

Y 

WSB 9 

16 21 34.69 

-26 12 26.9 

K5 

1.20 

12.88 

8.864 

5 

Y 

N 

WSB 12 

16 22 18.52 

-23 21 48.0 

K5 

1.20 

13.03 

8.109 

12 

Y 

N 

WSB 13 

16 22 22.31 

-25 52 19.9 

M4 

0.29 

14.56 

9.438 

5 

N 

N 


Note.— R magnitudes taken from NOMAD ^Zacharias et£^J200^), K magnitudes taken from 2MASS dSkrut^ie et aTill|2006[). For 
targets whose most recent classification in the literature proviaed a range of possible spectral types, the mid-range spectral tyjDe was 
used. 

Table[^is published in its entirety in the electronic edition of The Astrophysical Journal. A portion is shown here for guidance regarding 
its form and content. 
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by 

Simon et al. 

(1995 

I, 

Koresko (2002), Bar- 

Sony et al. 

(12003), 

Ratzka et al. 

(2005) and Cieza 

et al. (2010). The 0.9arcsec binary [MMG98J RX 


J1622.7-2325a was also added. This companion 
was found by Prato (2007), from AO images taken 
in conjunction with their spectroscopic observa¬ 
tions. Both components were spectrally classified, 
and so we adopt a contrast ratio that preserves this 
classification. In total, we found 40 companions to 
members of our target list (8 in triple systems and 
32 in binary systems). 


Due to difficulties with observing wide equal bi¬ 
nary systems the majority of known multiple sys¬ 
tems were not reobserved, which will introduce 
a bias against the detection of higher order sys¬ 
tems. The exceptions to this were the stars ROXs 
12, ROXs 42B, ROXs 47A, EM* SR 20, EM* SR 
21, EM* SR 24 S-N, GSS 31 and [MMG98] RX 
J1622.7-2325a. The secondaries in these systems 
were either wide enough to allow the individual 
components to be studied separately (separation 
> 4”), close enough to have little effect on the 
performance of the AO system (separation < 1”), 
or faint enough to allow the AO system to lock on 
the primary (contrast ratio < 0.5). The EM* SR 
24 S-N and ROXs 43A-B components were bright 
enough to make our target list individually and 
are included separately. 


Two of our targets host known wide compan¬ 
ions with masses close to or below the deuterium 
burning limit (ROXs 42B and ROXs 12). Both 
have been confirmed as comoving with their host 
star ( Kraus et al.||2014 Currie et al.||2014 |. 

The star DoAr 21 was also found to be a bi¬ 
nary in VLBA measurements by |Loinard et al.| 
(2008). However, the small angular separation of 
the components (5 mas) put it outside our region 
of interest and so it was reobserved and considered 
a single star in our analysis to avoid biasing our 
results. 


We have included the detection limits from the 
multiplicity surveys in the literature for both ob¬ 
served multiples and non-detections to provide in¬ 
formation on more widely separated companions. 
Combined, this data covers 80 of our targets. 

Two M dwarfs that failed our spectral type cut 
were observed: EM* SR 22 (M4.5) and SSTc2d 
J162224.4-245019 (M5). The stars [MMG98] RX 
J1625.2-2455b, [MMG98] RX J1624.8-2359 and 


WSB 74 were also observed. These were con¬ 
sidered as possible members by several surveys, 
but have not have their membership hrmly estab¬ 
lished. Finally, Elia 2-29 was observed, despite 
having no published spectral type. The results of 
these observations are reported here, but are not 
included in our statistical analysis or target list. 


2.1. Stellar and companion properties 


Stellar properties for pre-main sequence stars 
are difficult to estimate, due to a lack of observa¬ 
tional data to constrain stellar evolutionary mod¬ 
els. Masses can be uncertain by factors of 2 or 
more, due to unresolved multiplicity or intrinsic 
variability. For these reasons, any inferred proper¬ 
ties should be treated with caution. Relative prop¬ 
erties, such as the mass ratio q are less affected by 
these systematics and are much less uncertain. 

Masses for each target were estimated using the 


methods outlined in Kraus & Hillenbrand (20071. 


Since no single set of stellar models spans the en¬ 
tire range of spectral types found here, it was 
necessary to combine several independent mod¬ 
els. For high mass stars (spectral type F2 and 


earlier) the temperature scales of Schmidt-Kaler 


(19821 were used to estimate masses directly from 
the spectral type. For lower mass stars the ef¬ 
fective temperatures were estimated by combining 


the temperature scales of Schmidt-Kaler (1982) 


with the M dwarf temperature scales from Luh- 
man et al. (2003). The masses were then calcu- 


_^ ( |2003D . _ _ 

lated using the isochrones of |Chabrier et al. (2000) 
(COO), Baraffe et al. (1998) (B98) and Siess et al. 


(2000) (SOO) at IMyr, approximately the median 
age of Ophiuchus members. The COO isochrones 
were used for 0.001-0.1 M©, B98 was used for 0.1- 
0.5 Mq, a mass weighted average of B98 and SOO 
was used for the range 0.5-1.0Mq, and SOO was 
used for M > 1.0 q. 

Companion masses were calculated using the 
predicted absolute magnitudes from the isochrones 
and the measured contrast in the observed band. 

To determine which of our targets host disks, 
we use data from the “From Molecular Cores 
to Planet Forming Disks” Spitzer legacy survey 


(Evans et al. 2003), which covered most of Ophi¬ 


uchus. One product of this survey was the identifi¬ 
cation of YSO candidates based on the presence of 
infrared excesses across SED measurements cover- 
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ing 4.5-24 and we adopt these results. 
For those targets without a classification, 


we 


use the infrared excess measured by WISE (Wright 
et al.||2010|. We use an approach similar to Luh- 
man & Mamajek (2012), by comparing the Kg-WA 


colour as a function of spectral type for all targets. 
This metric easily separates the disk hosting pop¬ 
ulation. 


3. Observations and Data Analysis 
3.1. Observations 

The targets were observed over several years 
using the Keck NIRC2 instrument with its visible 
AO system, from 2008 through 2013. One dataset 
was also taken with the VLT NACO instrument 
and IR AO system, allowing targets with bright 
IR magnitudes to be observed. Both instruments 
have several aperture masks installed in a pupil 
wheel ( Tuthill et al.|[2006 Tuthill et al. 20101. A 
summary of the observations can be seen in Table 

H 

The observing strategy for our observations was 


similar to that used in Kraus et al. (2008). Targets 


were observed in groups of ~ 5 — 19 objects based 
on their brightness and location in the sky, with 
nearby calibrator stars regularly observed to es¬ 
timate instrument systematics. Calibrators were 
chosen from the 2MASS catalogue to have similar 
visible and infrared magnitudes to the targets in 
each group. In addition, targets that showed no 
detectable signal were used as calibrators for the 
remaining targets observed during the same night. 
In total 63 targets were observed, with total inte¬ 
gration times between 8—160 s based on the target 
brightness and chosen to give a similar number of 
total counts on the detector. 


3.2. Data Analysis and Detection Limits 

The aperture masking data was analysed with a 
data pipeline developed at the University of Syd¬ 
ney, one of two commonly used software packages 
for processing such data. A more thorough de¬ 
scription of this pipeline can be found in other 


work (e.g. Kraus et al. 2008), but a short summary 
follows. Images are background subtracted, flat 


^For futher details on the criteria used to identify disk host¬ 
ing stars in the Spitzer survey, see [Evans et al. ||2007||. 


fielded and windowed with a super-Gaussian func¬ 
tion with a FWHM of 500 mas, before complex 
visibilities are measured from the Fourier trans¬ 
form of the image. These raw visibilities are then 
turned into closure phases using a matched filter 
approach. Closure phase calibration is achieved 
by estimating the intrumental closure phases from 
the average of those measured on the calibrator 
stars, then subtracting this estimate from the tar¬ 
get closure phases. 


To determine the limits for companion detec¬ 
tion, a Monte-Carlo approach was taken. For each 
dataset a set of 10,000 random closure phases was 
drawn from a Gaussian distribution, with a width 
set by the uncertainty on each closure phase. A 
model fit was performed for each set of closure 
phases, yielding a list of 10,000 fake detections 
that were fit only to noise. To be considered bona- 
fide, a companion had to have a contrast ratio 
lower than a 99.9% {3.3a) detection limit, cal¬ 
culated as the contrast ratio for which 99.9% of 
the fake detections with a similar separation had 
a higher contrast. 


The detection limits calculated from each ob¬ 
servation are shown plotted in Figure showing 
the number of targets for which a companion with 
a given set of parameters could be detected. 

For the stars ROXs 47A and [MMG98] RX 
J1622.7-2325a, a different calibration scheme was 
used. ROXs 47A star is a known tertiary system 
from Barsony et al. (2003). The wider secondary 
companion was well within the instrument field of 
view, and appeared to have a significant influence 
on the measured closure phases of the primary. 
To account for this, the closure phases from the 
secondary component were measured and used to 
calibrate those from the primary. This resulted 
in a significant improvement in calibration, and a 
clean detection of the tertiary component. Sim¬ 
ilarly, [MMG98] RX J1622.7-2325a was a known 
0.9arcsec binary from Prato (2007), which placed 
both components well inside the instrument field 
of view. Inspection of the images showed that this 
system was in fact a triple, with a close pair and 
a brighter primary at a larger separation. Using 
the measured closure phases from the primary to 
calibrate those of the close pair produced a clean 
detection of the tertiary component. 
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Table 2 

Summary of aperture masking observations 


Date 

Telescope 

Aperture Mask 

Filter 

Targets Observed 

18/06/2008 

Keck 

9 hole 

Kp 

[MMG98] RX J1620.9-2352, [MMG98] RX J1621.4-2312, WSB 12, 
GSC 6794-537, 2MASS J16233234-2523485, DoAr 21, EM* SR 3, EM* 
SR 24 S, EM* SR 21, EM* SR 10, [MMG98] RX J1621.4-2332b, EM* 
SR 6, [MMG98] RX J1625.2-2455b 

31/05/2009 

Keck 

IS hole 

Hcont 

HIP 80126 

01/06/2009 

Keck 

9 hole 

CH4_short 

EM* SR 6, ROXs 47A 

05/04/2010 

Keck 

9 hole 

CH4_short 

EM* SR 24 S, ROXs 47A, Haro 1-16 

23/04/2011 

Keck 

9 hole 

Kp 

DoAr 25, EM* SR 8, WSB 12, DoAr 49, DoAr 24, WSB 46, WLY 2- 
10, DoAr 32, SSTc2d J162506.9-235050, WSB 63, SSTc2d J163355.6- 
244205, [MMG98] RX J1624.8-2359, WSB 74 

24/04/2011 

Keck 

9 hole 

Kp 

ROXs 3, DoAr 25, ROXs 12, GSS 31, WLY 2-10, DoAr 32, DoAr 33, 
[MMG98] RX J1623.8-2341a, ROXs 45D, SSTc2d J162506.9-235050, 
WSB 74, [MMG98] RX J1622.7-2325a 

04/06/2011 

Keck 

9 hole 

Lp 

DoAr 21, EM* SR 21, Haro 1-16 

05/06/2011 

Keck 

9 hole 

Lp 

DoAr 21, EM* SR 3, EM* SR 21, Haro 1-16 

23/06/2011 

Keck 

9 hole 

Kp 

EM* SR 24 S, ROXs 42B, EM* SR 6, EM* SR 20 

14/04/2012 

Keck 

9 hole 

CH4_short 

EM* SR 24 S, GSS 31, GSS 20, Haro 1-16, [MMG98] RX J1625.3-2402, 
WSB 40, [WMR2005] 1-21, EM* SR 22, EM* SR 20, GSS 35 

06/05/2012 

Keck 

9 hole 

Lp 

EM* SR 24 S, DoAr 25, EM* SR 21, GSS 31, DoAr 32, Haro 1-16 

07/07/2012 

Keck 

9 hole 

Kp 

WLY 2-10, [MMG98] RX J1622.6-2345, [MMG98] RX J1622.8-2333, 
EM* SR 6 

08/07/2012 

Keck 

9 hole 

Kp 

DoAr 25, GSS 31, DoAr 33, [MMG98] RX J1623.8-2341a, [MMG98] 
RX J1625.2-2455b 

09/03/2013 

VLT 

7 hole 

Ks 

SSTc2d J162224.4-245019, [E2011] 3-37, WLY 1-18, ROXs 39, GSS 
29, YLW 47, WSB 52, SSTc2d J162312.5-243641, WLY 1-13, GSS 26, 
[E2011] 6-62, SSTc2d J163603.9-242344, WSB 9, SSTc2d J162224.4- 
245019 

10/03/2013 

VLT 

7 hole 

Ks 

Elia 2-29, GSS 32, Elia 2-24, [E2011] 1-3, [E2011] 4-29, VSSG 
19, [MMG98] RX J1627.2-2404a , [MMG98] RX J1628.2-2405 , 
[WMR2005] 2-30, [MMG98] RX J1626.3-2407a , [E2011] 3-45, 

[WMR2005] 2-15, DoAr 13 

06/08/2013 

Keck 

9 hole 

Kp 

ROXs 42B 


Table 3 

Non Detections 

Primary 

Date 

Filter 

10-20 

20-40 

Am 

40-80 

80-160 

160-240 

240-320 

DoAr 21 

18/06/2008 

Kp 

2.08 

3.78 

4.65 

4.51 

4.32 

4.35 

DoAr 21 

18/06/2008 

Kp 

1.53 

3.41 

4.34 

4.03 

4.06 

4.06 

DoAr 21 

04/06/2011 

Lp 

0.0 

3.0 

4.84 

5.35 

5.11 

5.08 

DoAr 21 

05/06/2011 

Lp 

0.0 

3.18 

4.98 

5.49 

5.3 

5.26 

[MMG98] RX J1620.9-2352 

18/06/2008 

Kp 

2.15 

3.82 

4.69 

4.58 

4.38 

4.4 

[MMG98] RX J1620.9-2352 

18/06/2008 

Kp 

1.53 

3.41 

4.34 

4.05 

4.06 

4.06 

EM* SR 3 

18/06/2008 

Kp 

3.01 

4.53 

5.36 

5.28 

5.22 

5.28 

EM* SR 3 

05/06/2011 

Lp 

0.0 

2.81 

4.67 

5.17 

4.98 

4.95 

EM* SR 24 S 

18/06/2008 

Kp 

0.0 

0.82 

2.34 

1.92 

1.29 

0.92 

EM* SR 24 S 

18/06/2008 

Kp 

0.0 

1.4 

2.79 

2.45 

2.46 

2.46 


Note.— Angular separation ranges are given in mas. Tablef^is published in its entirety in the electronic edition of The Astrophysical 
Journal. A portion is shown here for guidance regarding its f^m and content. 
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Table 4 

Companions identified with aperture masking 


Primary 

Date 

Filter 

Separation 

(mas) 

P.A. 

(deg) 

Am 

q 

MJMp 

Mass 

Mg 

[MMG98] RX 

J1621.4-2332b^ 

18/06/2008 

Kp 

74.98 ± 0.42 

50.96 ± 0.30 

2.97 

± 0.03 

0.12 

0.10 

[MMG98] RX 

J1626.3-2407a^ 

10/03/2013 

Ks 

66.67 ± 0.45 

299.32 ± 0.45 

0.19 

± 0.02 

0.85 

0.33 

EM* SR 20 


23/06/2011 

Kp 

50.60 ± 0.30 

240.50 ± 0.30 

2.68 

± 0.02 

0.23 

0.83 

EM* SR 20 


14/04/2012 

CH4_short 

48.55 ± 0.31 

233.52 ± 0.36 

2.66 

± 0.03 

0.26 

0.92 

EM* SR 6“ 


18/06/2008 

Kp 

35.60 ± 0.20 

270.40 ± 0.15 

1.91 

± 0.01 

0.24 

0.72 

EM* SR 6“ 


01/06/2009 

GH4_short 

41.03 ± 0.22 

286.46 ± 0.33 

2.10 

± 0.02 

0.22 

0.66 

EM* SR 6“ 


23/06/2011 

Kp 

57.80 ± 0.20 

308.20 ± 0.22 

1.94 

± 0.02 

0.24 

0.70 

EM* SR 6“ 


07/07/2012 

Kp 

67.00 ± 1.20 

314.10 ± 0.95 

1.92 

± 0.07 

0.24 

0.71 

GSS 35 


14/04/2012 

Kp 

21.75 ± 0.39 

205.11 ± 0.74 

1.19 

± 0.19 

0.19 

1.45 

GSS 35 


09/03/2013 

Ks 

28.89 ± 3.47 

268.54 ± 2.13 

2.26 

± 0.43 

0.08 

0.60 

ROXs 39“ 


09/03/2013 

Ks 

46.17 ± 0.18 

191.10 ± 0.80 

0.20 

± 0.01 

0.86 

0.94 

ROXs 47A 


01/06/2009 

CH4_short 

51.50 ± 0.20 

99.20 ± 0.30 

0.19 

± 0.01 

0.90 

2.39 

ROXs 47A 


05/04/2010 

CH4_short 

43.38 ± 0.18 

109.06 ± 0.19 

0.22 

± 0.01 

0.88 

2.33 

[MMG98] RX 

J1622.7-2325a^ 

24/04/2011 

Kp 

128.16 ± 0.28 

280.53 ± 0.09 

0.14 

± 0.01 

0.90 

0.65 


^Companion reported for the first time in this work. 


4. Statistical framework 


To turn the measured detection limits and 
detections into estimates of the companion fre¬ 
quency distributions, we have employed a statis¬ 
tical framework similar to that of ICarson et al.l 


(2006), Lafreniere et al. (2007) and Evans et al. 


(2012). Further details may be found in these 


publications and the references therein, but a brief 
summary follows. 

If / is the fraction of stars with a companion in 
the range of masses [ruminj Wmax] and semi-major 
axes [amin^CLmax]^ & simple application of Bayes’ 
Theorem results in the following expression for the 
probability of / given the set of data {dj}: 




C{{d,}\f)P{f) 

C{{d,}\f)Pif)df 


( 1 ) 


Each dj is equal to 1 if a companion was detected 
and 0 otherwise. £{{dj}\f) is the likelihood of the 
data. We adopt a prior distribution on P{f) equal 
to the Jeffreys prior for a Bernoulli trial, P{f) = 
, ^ , but find no significant difference for any 

results when adopting a flat prior of P{f) = 1. 

If we let the probability of detecting a compan¬ 
ion in the given ranges of masses and semi-major 
axes, if it was present, be pj, then the probability 
of detecting a companion around any given star 
is fpj. Similarly, the probability of not detect¬ 
ing a companion around that star is 1 — fpj. The 


calculation of pj is described in section 4.1 


The likelihood of the data is then given by 
Equation 

( 2 ) 

We then calculate a credible interval for / us¬ 
ing the posterior P{f\{dj}), such that the inter¬ 
val between [fmimfmax] contains a fraction a of 
the total probability. We choose an equal tail dis¬ 
tribution, such that the probability contained in 
[0, fmin] and [fmax, 1] are equal. In the case where 
no companions are found in the stated interval, 
this is equivalent to solving Equation for an up¬ 
per bound on the companion fraction fmax- 


a = 



P{f\{dj})df 


( 3 ) 


If companions are found, the credible interval 
is found by solving the following two equations: 


P{f\{dm (4) 

^ P{f\{dm (5) 

fmax 

In this work we choose a conhdence level of 
a = 0.68 for our analysis, equivalent to l-cr limits, 
unless otherwise stated. 
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Fig. 3.— The companion detectability in differ¬ 
ent regions of the binary parameter space. Masses 
are calculated using the properties from Table 
and the approach described in Section [2.1[ Known 
binaries are marked with red triangles, while the 
new detections are shown with blue circles. The 
greyscale levels indicate the number of targets in 
our sample for which a companion with the given 
properties would have been detectable. 


4.1. Calculation of pj 


In order to calculate the probability of detect¬ 
ing existing companions in the intervals of mass 
and semi-major axes described above, we perform 
a Monte Carlo simulation. This involves generat¬ 
ing 10,000 companions with masses and angular 
separations according to known or proposed dis¬ 
tributions, and then comparing their properties 
to the calculated detection limits. The fraction 
of companions that fall above the mass detection 
limits is then used as an estimate of pj. 


When drawing masses for our simulated com¬ 
panions, we considered two distributions. First is 
the universal mass function of IMetchev &: Hillen-I 


brand (2009), proposed for companions to solar 


mass stars, given by 


dN 

dq 


oc q 


0.39 


( 6 ) 


We also consider a flat mass ratio distribu¬ 
tion, consistent with the observed data from |Kraus| 
et al. (2008) and Kraus et al. (2011). 


Our results indicate no significant difference in 
results between these distributions, and only the 
results of the flat mass ratio distribution are re¬ 
ported here. 

Angular separations are calculated by combin¬ 
ing information about semi-major axes, eccentrici¬ 
ties, orbital phases and inclinations. For the com¬ 
panion eccentricities, we follow the approach of 
Evans et aH (2012). In the absence of constraints 
on the eccentricity distribution of companions, we 
choose to draw them from a distribution of the 
form 

/(e) = 2e (7) 

which is arrived at based on theoretical consider¬ 
ations ( Ambartsumian|[l937 ). 

We convert the semi-major axis a into a pro¬ 
jected separation through multiplication by a pro¬ 
jection factor s, which is calculated from the ec¬ 
centricity distribution. Combining the above func¬ 
tion with the assumption that observing an or¬ 
bit from any direction is equally likely, [Brandeker 
et al. (2006) showed that the cumulative probabil¬ 


ity distribution for s is well approximated by the 
function 


E,(s) =0.5 


1 — cos I —s 


( 8 ) 
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We consider s in the interval [0, 2]. 

Finally, the semi-major axes of our simulated 
companions are drawn from a distribution of the 
form 


dN 

da 


( 9 ) 


This is consistent with the results of IKraus et al.l 


(2008, 2011). 


Since our approach involves explicitly assuming 
forms for the distributions of mass and separation, 
it can introduce systematic biases to the calculated 
companion fractions. However, since these apply 
only within each bin of mass and angular sepa¬ 
ration we can check our assumptions by compar¬ 
ing the observed distributions of these parameters 
across many bins to those used in our simulations. 


5. Stellar Multiplicity of Ophiuchus Mem¬ 
bers 

5.1. Detected Companions 

Our aperture masking observations were sensi¬ 
tive to stellar companions at separations between 
10-320 mas, while the imaging results from the lit¬ 
erature were sensitive to binary systems at up to 
Oarcsec. The detection limits for stars that were 
identified as single stars in our masking observa¬ 
tions are listed in Table [21 

Eight binary systems were identified from our 
observations, 5 for the first time, and are listed 
in Table The known close binary EM* SR 20 
was recovered in two epochs, with position angles 
that suggest at least one complete orbit since its 
earliest resolved observations in 1990. Both EM* 
SR 20 and the new multiple system EM* SR 6 
have enough epochs to formally allow an orbital 
fit. In addition, the star [MMG98] RX J1622.7- 
2325a was shown to be a triple system for the first 
time. 

While observing EM* SR 24 S, it was appar¬ 
ent from the AO guide camera that the contrast 
ratio between EM* SR 24 S and EM* SR 24 N 
has changed significantly over time. The S com¬ 
ponent provided much more flux than expected to 
the AO wavefront sensor and appeared noticeably 
brighter than the N component, suggesting that 
it has changed substantially since the epoch of its 
UCAC4 photometric measurements. Since this in¬ 
formation is not saved during NIRC2 observations, 


it is not possible to report a quantitative flux mea¬ 
surement. EM* SR 24 S is a known transition disk 
system, and similar variability is seen in many of 
these objects. 


Our observations failed to recover the close 


companion to ROXs 42B detected in Simon et al. 


(1995) and Ratzka et al. (2005) at separations of 
56 mas and 83 mas respectively. This companion 
should have been detectable at separations from 
10-320 mas, and by examining the raw images we 
can extend the outer limit to beyond 1 arcsec. The 
non-detection of this companion in our two obser¬ 
vation epochs separated by 2 years makes it un¬ 
likely that orbital motion may have carried it in¬ 
wards beyond our inner working angle. We have 
included the previously reported companion in our 
analysis, but note that further study is required to 
determine the nature of this system. 


None of the companions reported here have con¬ 
trasts or calculated masses consistent with brown 
dwarfs, despite 44 of our 63 masking observations 
reaching contrasts deep enough to detect such 
components. 


5.2. Dinarity Fraction 

By combining high resolution multiplicity data 
from the literature with our SAM observations, 
our survey covers binary and high mass substel- 
lar companions separated by 1.3-780 AU. Overall, 
we find 36 multiple systems, while 33 members 
of our target list are single stars across the full 
range of separations. The remaining 43 targets 
are missing either imaging or SAM observations, 
and so have incomplete information. If each ob¬ 
servation was able to detect all binaries across this 
separation range, these numbers suggest an overall 
binary fraction of 53 ± 6%, ignoring targets that 
were not observed with at least one of imaging or 
SAM. Similarly, when restricting the spatial scales 
to 1.3-41.6 AU to coincide with those explored by 
the SAM observations, we find 41 single stars and 
22 stars with one or more companions, yielding a 
binary fraction of 35 ± 6%. 

Using the Bayesian techniques described above 
we are able to include the completeness of cover¬ 
age of each observed target to better estimate the 
total companion frequency. Considering a range 
of separations between 1.3-780 AU and compan¬ 
ion masses greater than 0.08 Mq, we find a total 
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companion frequency of 43 ± 6 %. This value is 
between the companion frequency of found 


in Upper Scorpius (iKraus et al.|2008|) and 64^g^% 


in Taurus-Auriga (Kraus et al. 2011 


Previous surveys of Ophiuchus members have 
reported similar companion frequencies. |Ratzka| 
et al.j (20051 found that 29 ± 4% of targets were 
in multiple systems over the range 0.13-6.4 arcsec 
(17-830 AU), compared with our value of 32.3 ± 
5.3% over the same range. Our results reflect that 
of ||Ratzka et al. (2005), who found significant dif¬ 
ferences in stellar multiplicity between the Ophi¬ 
uchus and Taurus-Auriga star forming regions. 
Despite similarities between the two regions, the 
companion fraction of Ophiuchus is notably lower. 

Across the range of separations explored here, 
nine of our targets are part of systems with 3 com¬ 
ponents (EM* SR 13, EM* SR 24 N/S, WSB 18, 
WSB 38, SSTc2d J162944.3-244122, ROXs 42B, 
ROXs 43 A/B, ROXs 47A, [MMG98] RX J1622.7- 
2325a), but we find no evidence of higher order 
multiplicity across this separation range. 


5.3. The Mass Ratio Distribution 


A large survey of nearby field dwarfs by |Fis-| 
cher & Marcy (19921 found that the distribution 


of companion masses is consistent with a flat mass 
ratio distribution for ratios q > 0.4. In addi¬ 
tion, they found some evidence for a slight excess 
of equal mass binaries for M types, a conclusion 
echoed by Reid & Gizis ( |1997 ). This trend also ap¬ 
pears to hold for late M dwarfs and brown dwarfs 
( Glose et al.||2003 Bony et ar]|2003 1. In contrast, 
the results of Duquennoy & Mayor (1991) indi¬ 


cated a lack of equal mass binaries around F and 


G type stars. Similarly, Kraus et al. (2008) found a 


flat mass ratio distribution holds across the entire 
range of mass ratios for binary systems in Upper 
Scorpius, although with no evidence for a peak 
close to unity. 

In Figure we have shown the observed frac¬ 
tion of stars with companions with separations be¬ 
tween 1.3-780 AU as a function of their mass ratio. 
The observed distribution appears consistent with 
a flat mass ratio distribution. We find an excess 
of equal mass binaries at low significance, echoing 
the results of|Fischer & Marcy (1992) and Reid 


& Gizis (1997). However, we also find a similar 


excess of low mass ratio companions at low signif- 


0.35 

0.30 

0.25 

0.20 

0.15 

0.10 

0.05 

0.00 






0.2-10Mg 





1 ■ 


T , 1 


, 1 ■ 


- 1 - i 

1 1 1 - 

-J -1 

II III II 1 

1 It MM 1 II 

1 

J _ 1 III _ 

III. III! Ill 


gO.35 
0.30 
5 0.25 

c 0.20 

•§ 0.15 
0.10 
E 0.05 

S 0.00 


0.2-l.rLilU 


j_i _ I II I I 



Fig. 4.— The observed binary fraction of our tar¬ 
gets as a function of the mass ratio between the 
primary and secondary components, for the whole 
sample (top) as well as the low mass (< 1.1 Mq, 
middle) and high mass (> 1.1 M©, bottom) mem¬ 
bers. Gompanion separations between 1.3-780 AU 
were considered. The individual mass ratios are 
plotted as vertical lines on the x-axis. 


icance. 

A one-sided Kolmogorov-Smirnov test shows 
that the observed distribution is consistent with 
a flat mass ratio distribution, with test statistic 
D=0.15 and a p value of 0.21. 

Since the companion mass ratio distribution is 
expected to be dependent on the mass of the pri¬ 
mary star, we have shown the results of splitting 
our target list into a high mass and a low mass sub¬ 
sample. Using a cutoff of 1.1 Mq, approximately 
the median mass of our sample, the resulting mass 
ratio distributions are also shown in Figure]^ We 
find no significant difference between the shape of 
the two distributions, although the total compan¬ 
ion fraction appears lower for the high mass sam¬ 
ple. 
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5.4. The Separation Distribution 


Studies of binaries in the field have shown that 
the separation distribution (or equivalently, the 
period distribution) is approximately log-normal 
across a wide range of masses (e.g. Duquennoy 
fc Mayor 1991[ [Fischer &: Marc}^|1992y Raghavan 


et al. 


20101. In addition to this, the mean and 


standard deviation appear to be mass-dependent, 
with low mass stars having a lower mean separa¬ 


tion (Duquennoy & Mayor 1991 Close et al. 20031. 


However, binaries in young star forming regions 
may be better matched by a log-flat distribution 
with an outer cutoff that increases with the mass 
of the primary ( Kraus et al.||2008 20111. 

Figure shows the observed fraction of stars 
with companions with masses between 0.08-6 Mq 
as a function of the separation between the com¬ 
ponents for our sample. We have opted to plot 
the minimum observed separation between com¬ 
ponents (when multiple detection epochs or multi¬ 
ple companions are present), rather than the more 
physically meaningful semi-major axis of the or¬ 
bit. For the majority of our binary targets, not 
enough detection epochs or orbital motion has 
been observed to allow an orbital solution to be 
derived, and conversions between projected sepa¬ 
rations and semi-major axes rely on extrapolated 
properties of binaries with much smaller separa¬ 
tions. 


The observed distribution appears log-normal 
in shape, with a large standard deviation. Over 
the range of separations explored, there is lit¬ 
tle evidence to distinguish between the proposed 
log-normal and log-flat distributions. A one 
sided Kolmogorov-Smirnov test comparing the 
observed distribution with a log-flat distribution 
gives D=0.12 and p=0.46, while the same test on a 
log-normal distribution gives D=0.11 and p=0.66. 

Figure also shows the result of splitting our 
targets based on primary mass using a cutoff of 
1.1 Mq, approximately the median mass of the 
stars in our sample. The separation distribu¬ 
tion for the high mass members appears to follow 
a log-flat separation distribution, consistent with 
previous results from Upper Scorpius and Tau- 


rus ( 

Kraus et al.|2008 

Kraus & Hillenbrand|2009 

Kraus et al. 

2011 

). For the low mass subsample, 


the distribution appears closer to log-normal in 
shape, as we find no binary companions with sepa¬ 


rations between 200-780 AU and only 1 companion 
with a separation between 1.3-6 AU. This differ¬ 
ence between the companion distributions of high 
and low mass members is similar to that found in 


Taurus by Kraus et al. (20111. 


We can use Bayesian analysis techniques to 
place limits on possible distributions by comparing 
the observed separations to a log-normal distribu¬ 
tion with two parameters: the mean separation /i 
and the standard deviation a. This approach mir¬ 
rors that of similar analyses in the literature (e.g. 


Allen|2007 Kraus et al.|2008 Rizzuto et al.|2013| . 

The likelihood of the observed separations given 
particular values of and cr is given by: 




exp 

- (log ^ - log rif /2a^ 


rlog 

Jlog 

exp 

- (log fi - xf /2a^ 

dx 


( 10 ) 


where rmi„ and r„iax are the minimum and maxi¬ 
mum separations explored, set as 1.3-780 AU here. 


Using a uniform prior on /i and cr, we find that 
the most probable values for /r and cr for the full 
distribution are given by 50 and 1.0 respectively. 
However, the family of solutions approximating 
flat distributions (cr ^ 1, /c unbound) are within 
our uncertainties. The probability space for /x and 
cr is shown in Figure Restricting cr to be in the 
range [0.5, 2], we find that the mean separation is 
given by /X = 50)^20° AU, while the standard de¬ 
viation is cr = l.O^Q®. The resulting limits vary 
substantially based on the choice of prior distribu¬ 
tion. 


6. Substellar Companions 

Many of the detection limits of our aperture 
masking survey reach the brown dwarf mass 
regime (~13-80Mj), and combined with the de¬ 
tection limits estimated from the imaging data of 
previous surveys we can estimate the occurence 
rate of brown dwarf companions to our targets. 

Two of our targets have known widely sepa¬ 
rated companions that are close to the deuterium 
burning limit and have been confirmed as co¬ 
moving with their host stars. Our mass estimates 
for these companions places ROXs 42B b in the 
planetary regime (10 Mj) and ROXs 12 b in the 
brown dwarf regime (25Mj), close to the values 
of 10 ± 4Mjand 16 ± 4Mjcalculated by [Kraus 
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Fig. 5.— The observed binary fraction of our tar¬ 
gets as a function of the separation between the 
primary and secondary components for the whole 
sample (top) as well as the low mass (< 1.1 Mq, 
middle) and high mass (> 1.1 Mq, top) members. 
Secondary masses between 0.08-6.0 Mq, and sepa¬ 
rations between 1.3-780 AU were considered. The 
individual separations are plotted as vertical lines 
on the x-axis. 


(which incorporated the individual 
estimated ages of these systems). 

The companion to the star WSB 28 also has 
a mass consistent with a high mass brown dwarf 


(75 Mj) and is likely comoving. This companion 


was identified by 

Reipurth & Zinnecker 

(1993 

I and 

recovered by 

Ratzka et al. 

(2005 

)• 



We detect no additional brown dwarf com¬ 


panions to our targets. For companions with 
masses between 13-80 Mjand separations between 
1.3-780 AU, we find an observed frequency of 7^®% 
for our observed targets. 

For widely separated companions between 42- 
780 AU, we hnd the frequency to be 4^^% while for 
close companions with separations between 1.3- 
42 AU we obtain a 1-cr upper limit of 12%. 



Fig. 6.— The probability space for the mean sepa¬ 
ration (/i) and the standard deviation (u) assum¬ 
ing a log-normal distribution for the companion 
separations. Contours enclosing 25%, 50%, 75%, 
90%, 95% and 99% of the total probability are 
shown. The most likely values are /r = 50 AU and 
fT=1.0, but a log-flat distribution {a » 1) is not 
ruled out. 


7. The Effects of Multiplicity on Disk Evo¬ 
lution 

Combining our census on the binarity of our 
targets with information on which stars host disks 
provides a fantastic opportunity to study the ways 
in which multiplicity affects the evolution of disks 
around young stars. 


Since 17 of our targets were selected based on 


the infrared surveys of 

Luhman & Rieke ( 

1999) 

and Cieza et al. 

(2010 

, their inclusion may in- 


troduce a bias towards disk hosting stars with in¬ 
frared excesses and so we have not included them 
in the analysis presented in this section. 


We find that 44 of our targets host disks, while 
52 have no disk, and 1 has insufficient information 
to classify it. This gives an overall disk frequency 
of 46 ± 5%, including those stars with no multi¬ 
plicity information. However, this figure changes 
substantially depending on the binarity of the tar¬ 
get. For binary targets, the ratio is 50 ±8%, while 
for confirmed single stars it is 70 ± 9%. 

Stars with incomplete multiplicity information 
have a disk frequency of 24 ±7%. This low disk fre- 


13 






























































quency is likely due to the spatial distribution of 
those targets. Previous multiplicity surveys have 
targeted the central region of the Li688 cloud, 
which hosts the youngest stars in our sample. Evi¬ 
dence suggests that star formation has been ongo¬ 
ing in some regions of the Ophiuchus cloud com¬ 
plex for as long as 10 Myr (Martin et al. 19981, and 
the targets furthest from the L1688 central region 
show a correspondingly lower disk fraction. The 
spatial distribution of our targets and their disk 
hosting status are shown in Figure 



Fig. 8.— The fraction of observed binary stars 
that host disks as a function of the separation be¬ 
tween the host star and the binary companion. 


Figure shows the result of splitting the stars 
based on the separation of the binary compan¬ 
ions. From our survey we find that close bi¬ 
nary systems (<40 AU) are less likely to host disks 
than wider binaries and single stars. This implies 
that the presence of a close binary companion ei¬ 
ther speeds up the process of disk dispersal or in¬ 
hibits its formation. This conclusion echoes that 
of Kraus et al. (2012), who found that only ^1/3 
of stars with companions at separations smaller 
than 40 AU hosted disks at ages of only 1-2 Myr. 
At the 1-2 Myr age of Ophiuchus, we find a similar 
result. 


At the age of our sample, stars with a wide 
companion (>40 AU) have a similar disk frequency 
as single stars (69 ± 12% compared to 70 ± 9%). 
This implies that the presence of a wide compan¬ 
ion does not signicantly affect disk evolution at an 


age of ^1-2 Myr. 


The fraction of stars that host circumstellar 
disks at such young ages is of critical importance 
for the formation of giant planets. The core accre¬ 


tion mode of star formation (Pollack et al. 1996) 


requires several Myr to efficiently produce giant 
planets. Since only ~ 1/3 of close binary systems 
in star forming regions still possess disks at an 
age of only 1-2 Myr, this implies that giant planets 
around such systems would be rare. This supports 
early results from radial velocity and transit sur- 


veys ( 

Desidera & Barbieri||2007 

Roell et al.||2012 

Wang et al.||2014 

)■ 


A significant population of our targets have 
SED gaps indicative of a transition disk. 13 
transition disk candidates were identified by 
: WSB 12, SSTc2d J162506.9- 
DoAr 32, DoAr 33, WSB 63, 
SSTc2d J163355.6-244205, SSTc2d J162245.4- 
243124, WSB 38, EM* SR 9, SSTc2d J162944.3- 
244122, ROXs 42C and DoAr 21. In addition, 4 
more stars have been identified as transition disks: 


Cieza et al. (2010 


235050, DoAr 25, 


EM* SR 24 S, Haro 1-16 (D oAr 44), EM* SR 21 
and YLW 46 (Oph IRS 48) ( Andrews et al.|[^011 
Geers et aL]|2007 ). 


Of these 17 targets, 4 have close companions (< 
40 AU) and may be circumbinary rather than tran¬ 
sitional (DoAr 21, ROXs 42C, SSTc2d J162944.3- 
244122, WSB 38). SSTc2d J162944.3-244122 and 
WSB 38 are known triple systems with a close bi¬ 
nary and wider tertiary, and it is unclear which 
component hosts the disk. In addition, EM* SR 24 
S, SSTc2d J162245.4-243124, EM* SR 9 and EM* 
SR 21 have companions at wider separations that 
may not be the source of their disk gaps. However, 
the remaining 9 targets have no detected stellar or 
substellar companion. 


Many scenarios have been proposed for the ori¬ 
gin of the disk gaps in transition disks. The 
favoured scenario invokes dust clearing due to the 
presence of giant planets forming in the disk. How¬ 
ever, for the majority of the transition disk can¬ 
didates in our survey we find a lack of stellar or 
brown dwarf companions capable of causing such 
gaps. This lends weight to the giant planet forma¬ 
tion hypothesis, showing that most of these disks 
are not circumbinary in nature. 
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Fig. 7.— The spatial distribution of stars in our survey, labelled according to whether they possess circum- 
stellar disks. Targets with no disk are marked with red stars, those with disks are labelled with green circles 
and targets with insufficient information to classify them are marked with blue diamonds. A high fraction of 
stars in the highest density regions have disks, while the disk fraction is lower in the West and North-West. 
Stars in these regions are expected to be older than those in the core, and many have lost their disks. 
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8. Conclusion 

We have performed a multiplicity survey of the 
p Ophiuchus cloud complex, studying the occur¬ 
rence rates of stellar and sub-stellar companions. 
This information was combined with knowledge of 
circumstellar disks among our targets from Spitzer 
to investigate how these properties relate. Our re¬ 
sults point to several significant conclusions for the 
processes of star and giant planet formation. 


1. We have identified 36 multiple systems, in¬ 
cluding 9 systems with 3 components. Of 
these, 5 stellar companions were resolved for 
the first time. This gives a binarity frac¬ 
tion of 43±6%, intermediate between the re¬ 
sults of similar star forming regions in Up¬ 
per Scorpius ( Kraus et al.|2008 ) and Tauris- 
Auriga ( Kraus et al.||20H ). 


2. The observed distribution of companion 
masses is consistent with a flat mass ratio 
distribution when considering components 
with separations between 1.3-780 AU, while 
the separation distribution is consistent with 
a log-normal or log-flat distribution over the 
same range when considering companions 
with masses between 0.08-6.0 Mq. These 
results agree with previous surveys of star 
forming regions and the field. 


3. The fraction of disk hosting stars depends 
strongly on the existence of a close (<40 AU) 
stellar-mass companion. Only ~l/3 of stars 
in close binary systems still host circum¬ 
stellar disks after only 1-2 Myr, suggesting 
that the presence of such companions signif¬ 
icantly speeds up the process of disk disper¬ 
sal or inhibits their formation. The lack of 
disks around close binary systems suggests 
that giant planets formed by the slow pro¬ 
cess of core accretion should be rare around 
such systems, a prediction supported by the 
results of early RV and transit surveys. 


4. In contrast, a more widely separated com¬ 
panion appears to have no effect on the pres¬ 
ence of a disk at 1-2 Myr. We find that 
^ 70% of both single stars and wide binaries 
retain their circumstellar disks at the age of 
Ophiuchus members. 


The authors wish to recognize and acknowl¬ 
edge the very significant cultural role and rever¬ 
ence that the summit of Mauna Kea has always 
had within the indigenous Hawaiian community. 
We are most fortunate to have the opportunity to 
conduct observations from this mountain. 
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and the Millennium Science Initiative (Chilean 
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